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Abstract. The observed microlensing events towards the 
LMC do not have yet a coherent explanation. If they are 
due to Galactic halo objects, the nature of these objects is 
puzzling — half the halo in dark 0.5 M Q objects. On the 
other hand, traditional models of the LMC predict a self- 
lensing optical depth about an order of magnitude too low, 
although characteristics of some of the observed events fa- 
vor a self-lensing explanation. We present here two mod- 
els of the LMC taking into account the correlation between 
the mass of the stars and their velocity dispersion: a thin 
Mestel disk, and an ellipsoidal model. Both yield opti- 
cal depths, event rates, and event duration distributions 
compatible with the observations. The grounds for such 
models are discussed, as well as their observational conse- 
quences. 

Key words: : Galaxy: halo, kinematics and dynamics, 
stellar content - Cosmology: dark matter, gravitational 
lensing 



1. Introduction 

One of the greatest uncertainties in the interpretation of 
the microlensing events detected towards the Large and 
the Small Magellanic Clouds by the MACHO (Alcock et 
al. 1997a, |Alcock et al. 1997b ) and the EROS experiments 
(Renault et al. 1997, Palanque-Delabrouille et al. 1998) 
resides in the determination of the location and conse- 
quently of the nature of the lenses. 

The MACHO results can be converted into an optical 
depth of 2.9^0 g x 10~ 7 , while the two candidate events 
of EROS would account for ~ 0.8 x 10~ 7 . Combining both 
results leads to an average optical depth of 2.1+q | x 10~ 7 



(Bennett 199S) 



While this observed optical depth towards the LMC 
is too small to account for a full dark halo of compact 



objects (expected r = 5 x 10~ 7 ), it is much too large to be 
explained by the contribution from known populations of 
low-mass stars in the disks of the Milky Way and the LMC: 
based on star counts, the contribution from the disk of ou r 
Galaxy is found to be tmw ~ 8x 10~ 9 ( Gould et al. 1997 ), 
and the global rate from known stars in both galaxies has 
been estimated by several authors to be < 20% (Wu 1994, 
[5ahu 1994 , ie Rujula et al. 1995) of that observed. 

Moreover, the location of a few lenses has been de- 
termined, either through the non-detection of parallax ef- 
fect for 97-SMC-l ( [Palanquc-Dclabrouillc et al. 199S| ), or 
through the observation of the caustic crossing in the 
case of binary lens e s, for 98-SMC-l (A fonso et al. 1998, 
Albrow et al. 1999|, [Alcock et al. 1999|, Rhie et al. 1999, 



Udalski et al. 1999) or MACHO-LMC-9 (Bennett et al. 
1996). The case for the latter is still unclear however, since 
the measured lens velocity, incompatible with a Galactic 
halo lens unless the source happens to be a binary star 
with the appropriate separation, seems too low for a LMC 
disk star. 

It has been pointed out by several authors (Kerins & 
Evans 1999, [Pi Stefano 1999] ) that the fact that the ob- 
served binary lenses are all located in the Clouds might 
indicate that all observed events arise from Clouds deflec- 
tors. Indeed, if the two observed binary lenses with caustic 
crossings are in the Clouds, then the expected number of 
self-lensing caustic crossing events has to be greater than 
.53 at 90% C.L. A very conservative estimate of the frac- 



tion of caustic crossing events (Mao & Paczihsky 1991 
pi Stefano 199S ) implies that we expect more than 5 back- 
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ground self-lensing events, which is far from negligible as 
compared to the total number of observed events. 

The case for self-lensing in the LMC and SMC is thus 
still open, despite the "proof" (Gould 1995) that its con- 
tribution to the optical depth is too low. As shown in 
(Palanque-Delabrouille et al 1998), the SMC could be thick 
enough to yield an expected optical depth as large as 
~ 2 x 10 -7 , compatible with the current observations to- 
wards this line-of-sight. A more thorough exploration of 
possible LMC models therefore appears worthwhile. 
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Gould (1995) noticed that the optical depth associated 
to the self-lensing of a stellar disk, a good approximation 
for the bulk of the lmc matter, depends only on the ver- 
tical velocity dispersion aw of the stars: 



molecular clouds. The diffusion of stellar orbits caused by 
such gravitational inhomogcncitics results in older popu- 
lations being associated with larger velocity dispersions. 
Stars diffuse in velocity space. For the Milky Way, Wielen 
obtained a global diffusion coefficient 



= 2 



^sec 2 i 



(1) 



C„ 



6. x 1(T 7 km 2 s" 2 



where i is the inclination of the disk. Measurements of 
stellar velocities in the LMC are difficult because of its dis- 
tance. Only the brightest stars are resolved. Gould based 
his analysis on the CH objects studied by Cowley & Hartwick 
(1991). These carbon stars belong to the asymptotic gi- 
ant branch (AGB) of the HR diagram. With an apparent 
magnitude mboi ~ 14 — 15 and assuming a distance mod- 
ulus of 18.6, their bolometric luminosities are in the range 
2 — 5 x 10 3 Lq. A velocity dispersion of ~ 20 km/s led 
Gould to the conclusion that the optical depth for lmc 
self-lensing could not exceed 1.1 x 1CP 8 , well below the 
measured value of 1 — 3 x 10~ 7 . However, carbon stars are 
a peculiar population. They have left the main sequence 
(ms) on which they have spent at most ~ 10 — 12 Gyr, the 
age of the LMC. The average lifetime of a star on the MS 
decreases with the initial stellar mass m as 



(4) 



'1.2 x 10 10 yr) 



(m [M ])- 



(2) 



so that the oldest carbon stars observed today had an ini- 
tial mass larger than ~ 1 M@. Assuming a Salpeter IMF 
for the LMC (Geha et al. 1998) and adopting the star for- 
mation history derived by the same authors lead to the 
conclusion that approximatively 75 % of the sample ana- 
lyzed by Cowley & Hartwick (1991) has formed less than 2 
Gyr ago, during a recent stellar burst. The velocity disper- 
sion of 20 km/s refers therefore to young objects whereas 
the bulk of the LMC population is not dominated by young 
stars (Geha et al. 1998) and could be associated to much 
larger velocity dispersions. As a matter of fact, Hughes 
et al. (1991) have measured the radial velocities of long 
period variables (LPV) in the LMC and showed that while 
young and intermediate LPV's had a dispersion velocity 
in the range between 12 and 25 km/s, the old LPV's had 
a significantly larger dispersion velocity ranging from 31 
to 45 km/s. 

The increase of the velocity dispersion with the age 
of stars is a well-known fact. Bienayme et al. (1987) have 
modeled stellar populations in the Milky Way disk and 
showed that while aw ~ 6 km/s at birth, it increases up 
to 25 km/s for an age of 7 — 10 Gyr. That is why the 
characteristic thickness for O and B stars which formed 
very recently is ~ 200 pc whereas it reaches ~ 700 pc 
for solar type objects. Wielen (1977) has shown that the 
velocity dispersion a and the age r of a stellar population 
are related by 



C„ 



(3) 



where ao is the initial velocity dispersion. As time goes on, 
disk stars undergo dynamical shocks with, for instance, 



whereas he got a value 6 times smaller for vertical motions. 

We estimate therefore that the quoted value aw ~ 
20 km/s is only representative of the youngest stars and 
that the bulk of the LMC populations are older, hence with 
larger velocity dispersions. The favored star formation his- 
tory (model e of Geha et al. 1998) has a star formation 
rate which remains constant for 10 Gyr and then increases 
by a factor of three for the past 2 Gyr. The average stellar 
age is ~ 5 Gyr, quite similar to the age of the Milky Way 
disk populations. For the latter, aw increases from 6 to 25 
km/s as stars grow old. Scaling these results to the LMC 
implies that the oldest populations should reach there a 
dispersion velocity aw ~ 80 km/s. The resulting optical 
depth would be therefore an order of magnitude larger 
than Gould's estimate, in fairly good agreement now with 
the observations. 

At this point, a detailed investigation of the effect of an 
increased velocity dispersion on the optical depth for LMC 
self-lensing is mandatory. To start with, we have refined 
in section 2 Gould's description by introducing several 
isothermal stellar components. Depending on the stellar 
species, the velocity dispersion aw varies from 20 up to, at 
most, 80 km/s. This corresponds to a diffusion coefficient 
C a lying in the range from to 5 x 10 -7 km 2 s~ 2 yr -1 . The 
gravitational potential satisfies a Poisson equation. From 
the determination of the vertical distribution of stars, the 
self-lensing optical depth and the event rate are calculated, 
along with the distribution of event durations. For our ex- 
treme model, the optical depth is ~ 0.7 — 1 x 10~ 7 in the 
inner 3°. As a first approach, the variation of the surface 
mass density S with radius R is assumed to follow a Mes- 
tel profile (see section ||) . We have also investigated the 
tidal effects that our Galaxy generates on the stellar popu- 
lations of the LMC. Notice that in this thin-disk model, the 
vertical and radial directions are somewhat disconnected. 
At radius R, the stellar distribution is derived as if the 
LMC disk looked infinite. We have therefore considered in 
section 3 the case of a multi-component model where each 
stellar population is distributed according to an oblate 
spheroid, with significant flattening along the disk. The ve- 
locity dispersion aw and the flattening depend on the stel- 
lar species. In the extreme model, a w ax reaches 45 km/s, 
hence a diffusion coefficient C CT ~ 1.4 x 10~ 7 km 2 s~ 2 yr" 1 , 
not too far from the Milky Way value. Finally, in section 4, 
we show that adopting a larger value for the lmc rotation 
velocity results into an enhanced self-lensing rate. We sug- 
gest further strategies to test our models and we conclude 
in section [|. 
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2. The Large Magellanic Cloud as a Thin Mestel 
Disk 

2.1. Description of the model 

The bulk of the mass of the LMC lies in a nearly face-on 
disk, with an inclination of 27°. The Cloud exhibit a bar 
structure emitting about 10 to 15 % of the total luminos- 
ity, but which is not included in the models presented here. 
This disk extends at least to 8° from the center but may 
conceivably be larger. The mass in the inner ~ 4 — 5 kpc 
(1° = 0.92 kpc for a distance modulus of 18.6) is estimated 
by Hughes et al. (1991) to be < 6.2±1.5xl0 9 M from the 
dispersion of the radial velocities of Long Period Variables 
(LPV) stars. The virial mass is smaller: 4.4 x 10 9 M Q . 

Measurements of the velocity dispersion of planetary 
nebulae and star clusters up to ~ 8° radius as well as HI 
observations indicate that the LMC disk rotates. The ro- 
tation curve appears flat out to about 8° radius, and may 
remain constant out to 15°. As the rotation curve of the 
LMC does not show signs of a Keplerian falloff, the mass 
in the inner 8° radius is inferred (Schommer et al. 1992) 
to be < 1.5 ± 2.5 x 10 10 M Q . Because the lmc disk is 
nearly face-on, the rotation velocity is hard to determine 
accurately. It lies in the range between 50 and 80 km/s. 
From now on, we adopt the median value of Vc = 65 
km/s, slightly smaller than the value of 77 km/s favored 
by Schommer et al. (1992). Our LMC mass is therefore on 
the low side. Note that our approach is meant to be conser- 
vative. Should the LMC be more massive, our conclusions 
would be strengthened. As shown in section 4, a larger 
rotation velocity would actually result into an enhanced 
rate for gravitational self-lensing. 

The surface mass density at radius R follows a Mestel 
profile 



E(i?) = 156[M pc" 2 ] 



V c 



65 km s 



This implies a LMC mass of 



M(R) = 4.9 x 1O 9 M 



V c 



65 kms 



1 kpc 
R 



R 
5 kpc 



(5) 



(G) 



Increasing the rotation velocity from 65 to 80 km/s results 
in an increase of the mass in the inner 5 kpc from 4.9 to 
7.4 x 10 9 Mq. Both values are in fair agreement with the 
value 4.7 - 7.7 x 10 9 M© of Hughes et al. (1991). 

We have modeled the disk to contain several stellar 
populations. We assumed that a Salpeter distribution holds 
for the surface mass densities: 



requirement that the global surface mass density S(i?) is 
recovered once the contributions of the various stellar 
components are added together. The luminosity of a MS 
star is related to its mass through L[Lq] = (m[M ])' 3 , 
with (3 w 4 (Henry & McCarthy 1993). The stars of the 
LMC that are monitored by the microlensing experiments 
EROS and MACHO have an absolute magnitude smaller 
than My ~ 1- This corresponds to stellar masses ~ 2.2 — 
2.4 M Q . We have therefore extended our stellar distribu- 
tion up to 77i max = 2Mq. In our simple model, this mas- 
sive component mimics the brightest members of the LMC 
populations, such as the CH stars. It also stands for the 
sources which may be potentially amplified by other LMC 
members. Notice that our results are not sensitive to the 
exact value of m max . The lightest stars of our model have a 
mass m m i n = 0.1 Mq, They correspond to the faintest ob- 
jects on the MS, the red-dwarves. The mass-to-light ratio 
of our stellar distribution is readily derived 



dm 



oc TO 



(7) 



where a — 1.35. Above 0.5 Mq, this is consistent with 
HST stellar counts in the LMC as discussed by Geha et 
al. (1998). The normalization factor is obtained from the 



M 
T 



0-a 
a- 1 



(8) 



With the values a = 1.35 and j3 = 4, this implies M/L = 
1.75 in solar units. At a radius of 1 kpc, the surface bright- 
ness is ~ 90 Lq pc~ 2 . Relation || gives a surface mass den- 
sity of 156M0pc -2 , hence a mass-to-light ratio of 1.71, 
in good agreement with our model (the effect of extinc- 
tion is neglected). Note that in the case of the oblate 
spheroids of section 3, the surface mass density E(l kpc) ~ 
220M Q pc -2 and the mass-to-light ratio becomes 2.4. Those 
values may be compared to the mass-to-light ratio of our 
galactic disk where M/L ~ 2. As regards the LMC, the 
general trend is an increase of M/L with the radius R. At 
4 and 5 kpc, the thin Mestel disk model respectively gives 
M/L — 4.3 and 5.3 whereas, for the flattened spheroids, 
we get 5 and 6. This is in agreement with Schommer et 
al. (1992) who conclude that the integrated mass-to-light 
ratio of the LMC is ~ 10. That increase of M/L with R 
might be due to an evolution of the stellar composition 
along the disk. The farther away from the center a stellar 
field lies, the fewer bright stars it contains. This should 
lead to different HR diagrams at different locations in the 
LMC disk. For simplicity, we have not dealt here with such 
refinements and we have kept our stellar distributions in- 
dependent of the position R. We leave to section 4 the 
discussion of the influence of M/L on our results. 

The vertical velocity dispersion aw depends on the 
stellar population. The only relevant dynamical parame- 
ter is actually the relationship between the surface mass 
density and aw- Stars of different masses may well have 
the same velocity dispersions. If so, their vertical distribu- 
tions are similar, with identical disk thickness. Relation |^ 
implies that stars with the same age should have the same 
velocity dispersion. As the LMC formed ~ 12 Gyr ago, stars 
whose mass to < 1 M have not yet left the MS. They have 
the same average age of 6 Gyr. Their velocity dispersion 
had ample time to diffuse and increase. That is why our 
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Table 1. Model A corresponds to the radius R = 1 kpc, 
a rotation velocity of 65 km/s and a w ax = 80 km/s. 



m[M @ ] 


aw [km / s] 


10 4 x p c [M Q pc- 6 ] 


£[M pc^] 


0.1 - 1 


80 


226 


134.4 


1.1 


68.5 


6.73 


3.3 


1.2 


59.2 


7.13 


2.9 


1.3 


51.6 


7.52 


2.6 


1.4 


45.2 


7.91 


2.4 


1.5 


39.7 


8.33 


2.2 


1.6 


34.9 


8.78 


2.0 


1.7 


30.6 


9.31 


1.8 


1.8 


26.8 


9.93 


1.7 


1.9 


23.3 


10.7 


1.6 


2 


20 


11.7 


1.5 



first stellar bin comprises stars with masses between 0.1 
Mq (red-dwarves) and 1 M© (solar type objects). Their 



velocity dispersion cr w &x is the largest of all species. Be- 
tween 1 and 2 M© (our most massive bin), the dispersion 
aw is obtained from the diffusion equation (0) and from 
noticing that the MS lifetime decreases as m . This leads 
to the relation 



'W 



A + 



(9) 



The coefficients A and B are determined by requiring that 
the dispersion velocity increases from <jq = 20 km/s for 
the brightest and heaviest objects up to a w ax for the light 
stars. Table 1 features the velocity dispersions of the var- 
ious stellar components in the case of model A where the 
radius R — 1 kpc and <r^ ax = 80 km/s. In the sequel, we 
will assume an horizontal velocity dispersion equal to the 
vertical velocity dispersion. 

Each model is specified by the radius R that has been 
varied from 1 to 5 kpc. The regions surveyed by EROS and 
MACHO are actually in the inner 5° from the LMC center. 
The value of a w &K has been varied from 20 to 80 km/s and 
the rotation velocity has been set equal to Vc = 65 km/s. 
Once a model is specified, the distribution of surface mass 
densities and vertical velocity dispersions ensues from the 
arguments discussed above. Each stellar component i is 
assumed to be isothermal, with vertical velocity dispersion 
aw,i- It is distributed according to 



Pi 0) 



Pi,c exp • 



'W,i 



(10) 



We have solved the Poisson equation along the vertical 
direction z 



dz 2 



4irGp(z) , 



where the total mass density p(z) comprises the various 
stellar species. The gravitational potential Q(z) and its 



derivative <I>'(z) are set equal to at the origin. The Pois- 
son equation is integrated up to z max = ±10 kpc. Stars 
located beyond that limit cannot really be considered as 
members of the lmc system. For each stellar population 
i, we compute the surface mass density for \z\ < 10 kpc 
and modify its central mass density p; iC in such a way that, 
after a dozen iterations, the proper distribution of surface 
mass densities is recovered. The convergence is found to 
be excellent. 

Should the lmc disk contain a single stellar species 
and extend vertically to infinity, we would have obtained 
a surface mass density of 



£ = 2v / 2 / o c a , 
where the typical thickness a is defined as 



V^Gpc 



(12) 



(13) 



Assuming a Mestcl profile to accommodate the flat rota- 
tion curve of the LMC would have led to the conclusion 
that the thickness at radius R is 

, 2 ' \ -2 

1.34 x 10 2 pc 



20 km/s J V 65 km/s 



R 



1 kpc 



(14) 



This simple picture is complicated by the presence of sev- 
eral stellar populations with different velocity dispersions 
and also by the presence of a cut-off at ± 10 kpc. As is 
clear from table [l| the low-mass stars with m = 0.1 — 1 Mq 
dominate the disk dynamics. The optical depth r LMC for 
self-lensing ensues from relation ([!]) where a w &x stands for 
the velocity dispersion. However, Gould's limit is only re- 
covered in the LMC regions where the disk thickness a is 
much smaller than the vertical boundary z max = ±10 kpc. 
This is the case for R = 1 kpc and for <r w a * — 20 km/s 
where a = 130 pc. When the velocity dispersion of low- 
mass stars is increased up to 80 km/s or when the radius R 
is large, the LMC disk is no longer thin with respect to the 
boundary z max . The actual value of r LMC is smaller than 
what would naively be inferred from a mere scaling from 
an infinite disk. This will be further discussed in section 4. 



2.2. Microlensing parameters 

For each of the models, one can compute the total optical 
depth t and the event rate T. The contribution of sources 
of mass m s located at z s and deflectors of mass md to 
these microlensing parameters is given by: 



(H) ST(m d ,z s ) 



dz. 



p d (z d cos 9, m d ) 



m d 



irR E (m d ,z d ,z s ) 2 (15) 
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5T(m d ,m s ,z s ) = 

Pd(zd cos 6>,m d ) 



dzd ■ 



m d 



2R E {m d , Zd, Z 8 ) <|«rol|> (16) 



where is the inclination angle of the LMC (taken equal 
to 27 degrees), z is measured along the line-of-sight, and 
p d is the mass density of the deflectors. 

The velocities v s and v d of the objects are supposed to 
be Gaussian-distributed with a mass-dependent velocity 
dispersion, as explained in section 2.1. The term < |u re i| > 
in equation ( |T6| ) represents the average of the relative ve- 
locity between the deflector and the source over all possi- 
ble velocities of these two objects: 

<M >=///.(«., m s )f d (v d ,m d )\v s -v d \dv s dv d 

(17) 

where f s (resp. f d ) is the source (resp. deflector) velocity 
distribution function — projection effects are neglected 
since z <C A,mc- 

Let p s be the mass density of the sources. The mass 
density at a given coordinate only depends on the mass 
of the object. Integrating equations (|l5|) and ( pf ) over all 
source positions yields: 



T(m d ,m s ) 



T(m d ,m s ) 



dz 



p s (z s cos 0.rn s ) 



5T(m d ,z s ) 



dz< 



p s (z s cos 8,m s ) 



dz 



p s (z s cos &,m s ) 



6T(m d , z s 



r. 



dz, 



p s (z& cos 0,rn s ) 



(18) 



(19) 



The values obtained with the equations (|18| ) and ( |19[) 
must then be summed on all the deflector species, and av- 
eraged on all the source species. Only the uppermost mass 
bin will be considered as possible sources (i.e. m s E [1.9 — 
—2.0] Mq) in agreement with the mass of the faintest stars 
that are reconstructed in microlensing surveys, whereas all 
the mass bins are of course taken into account when sum- 
ming for all possible deflectors. The global optical depth 
and event rate due to the entire stellar population are 
therefore given by: 



max )) (20) 

»d=Wx 

max )) (21) 

The first mass bin has in practice to be treated specially, 
since it is not thin. The deflector masses in this bin follow 
the Salpeter law given in Eq. (0). 

These results depend on the maximum dispersion ve- 
locity of the stars (i.e. on the dispersion velocity of the 
lightest population) , as illustrated in figures |l| and ||. The 
configuration where all the stellar populations have the 



= Ikpc 



V M „(in kmls] 



Fig. 1. Evolution of the optical depth r with the max- 
imum velocity dispersion of the various populations for 
the disk model, at two different distances from the center 
of the LMC: 1 kpc and 5 kpc. 



- Ikpc 



r = 5kpc 



V MAX ( in km,S > 



Fig. 2. Evolution of the event rate T with the maximum 
velocity dispersion of the various populations for the disk 
model, at two different distances from the center of the 
LMC: 1 kpc and 5 kpc. 



same velocity dispersion (i.e. when w max = 20 km/s) yields 
the same results as the calculations done by Gould (Gould 
1995). The slower increase in the optical depth with the 
maximum velocity for a radial distance of 5 kpc as com- 
pared to that observed at a distance of 1 kpc is due to 
the 10 kpc vertical cut-off we have imposed on the stars 
to be considered as LMC objects (see comment at the end 
of section 2.1), since the disk gets wider on the edges. As 
expected, the optical depth and event rate decrease as one 
goes to larger and larger radii from the galactic center. For 
?j max = 80 km/s, there is a factor of two in optical depths 
and three in event rates between a distance of 1 and 5 kpc. 

An interesting comparison to the microlensing data 
can be done through the predicted distribution of event 
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durations, P(At) which can be calculated directly from 
the event rate: P(Af) oc dT/dAt. In a similar manner as 
what was done for the calculation of the optical depth and 
the event rate, we have: 



dr f A 

-j^(At,m d ,m s ,z s 



Pd(z d cos9,m d ) 
dz d x 

dvddv s f d (v d ,m d ) f s (v s ,m s ) x 
Re 



(22) 



At 2R E \v s 
\v s -v d \ J 

Integrating over the source positions yields: 



v d \ 



dAt 



(At,m d ,m s ) 



p s (z s cos 6,m s ) dT 
m s dAt 



■j& t {At,m d ,m s ,z s ) 



(23) 



z s cos u,m s 



and therefore 



dAl( At ^= -^(At : m d (i d ),m s (i max )) 



(24) 



Figures |3| and ^ show the distribution of the event dura- 
tion obtained in this multi-mass disk model of the LMC 
normalized to a maximum value of 1, for two different 
values of the maximum velocity dispersion and at two dif- 
ferent distances from the galactic center. 



r = 1 kpc 






v mx = 2< " m/s 




v MAX = 80km/S 







Duration At (in days) 

Fig. 3. Predicted distribution of event durations for the 
multi-mass disk model at a distance of 1 kpc from the 
center of the LMC. Plain curve is for a maximum veloc- 
ity dispersion of 80 km/s, dashed curve for a maximum 
velocity dispersion of 20 km/s. 



At 1 kpc from the center, the duration distribution is 
nearly independent of the maximum velocity dispersion, 



r = 5 kpc 




Duration At (in days) 

Fig. 4. Predicted distribution of event durations for the 
multi-mass disk model at a distance of 5 kpc from the 
center of the LMC. Plain curve is for a maximum veloc- 
ity dispersion of 80 km/s, dashed curve for a maximum 
velocity dispersion of 20 km/s. 



because the velocity cancels out in Re/v when z oc v 2 . 
At 5 kpc from the center, the truncation at z = 10 kpc 
prevents this cancellation. 



2.3. Tidal effects 

We have finally modeled the tidal effects which the Milky 
Way potentially induces on the optical depth r LMC . Be- 
cause the LMC rotates in the potential well of our Galaxy, 
the interplay between the gravitational attraction and the 
centrifugal force results into an effective tidal potential 



3 GMr 



2 D 3 



- cos 2 i z 1 



(25) 



in the vicinity of the disk, for small z. The associated 
repulsion tends to stretch the stellar populations apart 
from the disk as if they were no longer pressure supported. 
The Poisson equation now becomes 



d 2 $ 

dz 



33 = AnG {p{z) - pg a l} 



(26) 



where p ga i denotes the average density of our Galaxy as 
seen by the LMC at a distance of D hMC = 52 kpc. With 
a galactic mass of M ga i = 6 x 10 11 M@, we get p ga i ~ 
10 -3 Mq pc -3 . The tidal forces make the potential well of 
the lmc disk flatter. 

In figure |5| the optical depth r LMC is plotted as a func- 
tion of the velocity dispersion crj^ ax , for various radii R. 
The trend is an increase of t lmc when tidal effects are 
included. Near the LMC center, the disk is thin with a 
strong cohesion. Tidal effects are negligible. Far from the 
LMC center, the vertical scale length of the disk increases, 



Aubourg et al.: Microlensing Optical Depth of the lmc 



7 



translating into an enhanced sensitivity of the stellar pop- 
ulations to the Milky Way tide. The increase of the optical 
depth is noticeable. Finally, for R > 5 kpc, the disk scale 
length a compares with the cut-off z max = ±10 kpc as 
the dispersion velocity cr^ ax increases towards 80 km/s. In 
that case, low-mass stars tend to be uniformly distributed, 
even in the absence of tidal forces. 





with tide effects 






without tide effects 


r = l kpc// 






/ y.~'"t = 3 kpe 












t = 7 kpc 



10 20 30 40 SO 



70 80 90 
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Fig. 5. Optical depth r LMC as a function of the velocity 
dispersion <7^ ax , for various radii R, with and without the 
inclusion of tidal effects in the model. 




10 15 



Fig. 6. Truncation ellipsoids, each corresponding to a dif- 
ferent species. Outer ellipsoid corresponds to the 0.1 Mq 
- 1 Mq mass bin, and then inwards to the species from 
1.1 to 2 Mq, by 0.1 Mq increments. 



The surface mass density of species i is 
Sj(-R) = — ^— '- avctan 



2 

max 



R 2 



R 



R 



(29) 



3. Ellipsoidal models 

3.1. Description of the model 

In the previous disk model, the z and R dependences of 
the potential were decoupled, which is not true for the 
exact potential of the lmc. In this section, we therefore 
investigate another model, more realistic as far as the co- 
dependence on z and R is concerned. 

The lmc is described with the same mass bins as in 
the previous section. The density of each population is 
assumed to be of the form 

Pi (R,z) = ^-j- (27) 

R 2 



so that the iso-density surfaces are oblate ellipsoids with 
ellipticity e,. This distribution is truncated at the ellipsoid 
passing through the cut-off radius i? max = 15 kpc, beyond 
which the mass density is zero. Note that the cut-off ellip- 
soids are not the same for the different mass species. They 
are illustrated in figure ^. 

The corresponding rotation curve is flat, with a circu- 
lar velocity 



AttG 



arcsm e» 



(28) 



With a given set of parameters (Aj,ej), the local vertical 
velocity dispersion is defined as 



a W,i 



$(R,z)-<f>(R,0) 
log{p(R,0)/p(R,z)} 



(30) 



We define vw,i as the vertical velocity dispersions aver- 
aged along z, which are more closely related to observable 
quantities. 



vw,i{R) 



J-oo Pi{R, z)a Wti (R, z)dz 

I-oo Pi( R ' Z ) dz 



(31) 



We computed the parameters Aj and e, that realize 
the following properties: 

• The circular velocity is 65 km/s. 

• The surface density mass function is given by a Salpeter 
law (see section 2). As can be seen in (^9|), this condi- 
tion can be consistently fulfilled for every radius. 

• as in section 2, the average vertical velocity dispersions 
vw,i vary according to the mass m, as given in (j^), from 

^min tO ^max* 

Unlike the disk model, there is an upper limit to the 
vertical velocity dispersion we can impose. For a one- 
species model, this limit would be cr^ ax = V c /V2, nu- 
merically 46 km/s, given the chosen rotation velocity. This 
limit is slightly different for a multi-mass model, but still 
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present. Accordingly, we impose velocity dispersions rang- 
ing from v m in = 20km/s for the heavy stars to v max = 
45km/s for the light stars. This is model B-l. For com- 
parison, we also present the results for model B-0 in which 
all species have a velocity dispersion of 20 km/s. The pa- 
rameters for model B-l are shown in table @. 



Table 2. Parameters obtained after convergence for model 
B-l, i.e. with v max = 45 km/s. 



m[M Q ] 




10^ x AifA/opc- 1 ] 


WK,i[km/s] 


0.1 - 1 


0.438 


70.5 


45 


1.1 


0.663 


2.08 


39.5 


1.2 


0.768 


2.16 


35.2 


1.3 


0.830 


2.22 


31.8 


1.4 


0.869 


2.27 


29.1 


1.5 


0.895 


2.29 


26.8 


1.6 


0.913 


2.30 


24.9 


1.7 


0.927 


2.30 


23.4 


1.8 


0.937 


2.29 


22.1 


1.9 


0.945 


2.26 


20.9 


2 


0.951 


2.23 


20 



3.2. Microlensing parameters 



In a similar way as that explained in section 2.2, it is pos- 
sible to compute the optical depth and event rate for an 
ellipsoidal model of the lmc. They are shown in figures 
and H as a function of the distance to the center of the 
LMC. There is a strong decrease of both quantities with 




Radial distance (in kpc) 



Fig. 8. Evolution of the event rate T with the distance 
to the center of the lmc, for the ellipsoidal model and 
two values of the maximum velocity dispersion: 20 km/s 
(model B-0) and 45 km/s (model B-l). 



the radius, which can be used to test statistically the lo- 
cation of the lenses detected in microlensing experiments. 
Note that the effect is even more drastic in the ellipsoidal 
model than in the disk model of the LMC: in the latter, 
the increased thickness at large radius was compensating 
the decrease in surface density. 

The predicted distributions of event durations, normal- 
ized to a maximum value of 1 , are shown in figures |9| and 
|i"0| at two different distances from the center of the LMC. 
Since there is no cut-off in the ellipsoidal model, the du- 
ration distribution is nearly independent of the maximum 
velocity dispersion at all radii. 




Radial distance (in kpc) 



Fig. 7. Evolution of the optical depth r with the distance 
to the center of the LMC, for the ellipsoidal model and 
two values of the maximum velocity dispersion: 20 km/s 
(model B-0) and 45 km/s (model B-l). 
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Fig. 9. Predicted distribution of event durations for the 
multi-mass ellipsoidal model at a distance of 1 kpc from 
the center of the lmc. Plain curve is for model B-l, dashed 
curve for model B-0. 
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r = 5 kpc 
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Fig. 10. Predicted distribution of event durations for the 
multi-mass ellipsoidal model at a distance of 5 kpc from 
the center of the LMC. Plain curve is for model B-l, dashed 
curve for model B-0. 



4. Discussion and conclusion 

We have presented here two models for the LMC, both con- 
sistent with the known properties of the Cloud (rotation 
velocity, surface light density), allowing for a multi-mass 
stellar population, each associated with a specific veloc- 
ity dispersion. This correlation between mass and velocity 
dispersion is justified by the fact that lighter stars are on 
the average older and had the time to diffuse in velocity 
space. 

These models of the LMC predict an optical depth 
in agreement with that observed by the MACHO and the 
EROS collaborations: at a distance of 1 kpc from the LMC 
center, the disk model yields r LMC (disk) ~ 1 x 10~ 7 for 
a velocity dispersion of the lightest and therefore old- 
est objects of 80 km/s, and the ellipsoidal model yields 
t L mo (ellipsoids) ~ 1.3 x 1CP 7 for w max = 45 km/s. 

The optical depth which we have derived results from 
the interplay of various effects. In the case of the disk 
model, for instance, Gould's estimate is only recovered in 
the thin disk limit. Setting the various stellar velocities 
at the same value of 80 km/s and pushing the LMC at 
infinite distance from the Milky Way, we get t lmc = 1.80 x 
10~ 7 , in agreement with relation (Q). If now we assume 
the same distribution of stellar velocities as in Table § 
i.e. with velocities ranging from 20 to 80 km/s, the optical 
depth becomes 1.14xl0~ 7 . Because the sources have small 
velocity dispersions, they are fairly concentrated towards 
the disk, hence a decrease of r LMC by a factor of ~ log 2. 
The proeminent effect is actually the presence of the cut- 
off z max = ±10 kpc which we have enforced. At 3 kpc from 
the center, the optical depth decreases from 1.14 x 10 -7 
to 7.57 x 10 -8 . At small radii, the surface mass density 
is large and the disk thickness a is small as compared 
to Zjaax- There, the decrease of t lmc is less sensitive than 



towards the outer fringes of the LMC disk, where the cutoff 
suppression is quite noticeable. Taking furthermore into 
account the finite distance of 52 kpc to the LMC implies 
a minor decrease of the optical depth which becomes now 
r LMC (3 kpc) = 6.63 x 10~ 8 . Finally, tidal forces tend to 
stretch the disk, hence a final value of r LMC (3 kpc) = 7.15 x 
10~ 8 . 

If the rotation velocity of the LMC Vc is set equal 
to 80 km/s instead of 65, the main impact on our disk 
model is an increased surface mass density. As the disk 
becomes thinner, the suppression of the optical depth re- 
sulting from the above mentioned cut-off is less sensitive, 
hence an increased optical depth at large radii. At 3 kpc, 
the optical depth becomes 9.08 x 10~ 8 . 

In order to model various M/L ratios, we have sim- 
ply modified the relative contribution of our first stellar 
bin. That population corresponds to faint stars. Still with 
Vc — 80 km/s and at R = 3 kpc, an increase of the M/L 
ratio by a factor of 2 and 5 respectively leads to 9.43 and 
9.67 x 10~ 8 . The optical depth is thus not sensitive to the 
precise value of the M/L ratio, insofar as the disk dynam- 
ics is dominated by the faint and low-mass species. 

Ellipsoids have a lenticular shape. As the radius R in- 
creases, their surface mass density drops, and the system 
becomes thinner towards the edges. Increasing the LMC 
rotation velocity without changing the velocity dispersion 
of the various stellar species does not change much the 
optical depth. Two opposite effects are actually at stake: 
the first is that the mass densities are larger, the second 
that the stellar populations are flattened towards the disk. 
If now the dispersion velocities aw are scaled to the rota- 
tion velocity, in order to maintain the shape of the ellip- 
soids — which depends to first order on the ratio aw /Vc 
— increasing the rotation speed from 65 to 80 km/s re- 
sults into a significantly larger optical depth. If we set 
°wr x = 55 km/s and Vc — 80 km/s, the optical depth at 
1 and 3 kpc becomes respectively 2.19 and 1.26 x 10~ 7 . 
That set of values actually correspond to an LMC mass at 
5 kpc of 7 x 10 9 M Q (Hughes et al. 1991) that increases up 
to 1.4 x 10 10 Af Q at 10 kpc (Schommer et al. 1992). 

Reexamining the status of event MACHO-LMC-9 ap- 
pears worthwhile. As pointed out in (Bennett, 1996), its 
measured projected relative velocity was too low both for 
a Galactic halo event and for self-lensing with a standard 
lmc disk model. Figure [ll] shows the distribution of the 
expected projected relative velocity for the 80 km/s disk 
and the 45 km/s ellipsoidal models, along with the mea- 
sured velocity, for two values of the horizontal velocity 
dispersion, null or equal to the vertical one - the reality 
being probably inbetween. The measured value is compat- 
ible with an ellipsoidal model, but only compatible with a 
disk model if the horizontal velocity dispersion is not too 
high. 

Reconsidering the shape of the LMC and the velocity 
distributions of its stellar components thus appears to be 
quite valuable. Two observational projects could yield in- 
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disk model 
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